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Abstract 

The transition between the two phases of 4D Euclidean Dynamical Triangulation [ ] was long 
believed to be of second order until in 1996 first order behavior was found for sufficiently large 
systems [5,9]. However, one may wonder if this finding was affected by the numerical methods 
used: to control volume fluctuations, in both studies [5,9] an artificial harmonic potential was 
added to the action and in [9] measurements were taken after a fixed number of accepted instead 
of attempted moves which introduces an additional error. Finally the simulations suffer from 
strong critical slowing down which may have been underestimated. 

In the present work, we address the above weaknesses: we allow the volume to fluctuate freely 
within a fixed interval; we take measurements after a fixed number of attempted moves; and we 
overcome critical slowing down by using an optimized parallel tempering algorithm [12]. With 
these improved methods, on systems of size up to N 4 = 64k 4-simplices, we confirm that the 
phase transition is 1 st order. 

In addition, we discuss a local criterion to decide whether parts of a triangulation are in the elon¬ 
gated or crumpled state and describe a new correspondence between EDT and the balls in boxes 
model. The latter gives rise to a modified partition function with an additional, third coupling. 
Finally, we propose and motivate a class of modified path-integral measures that might remove 
the metastability of the Markov chain and turn the phase transition into 2 nd order. 
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1 Introduction 


Euclidean Dynamical Triangulation (EDT) in four dimensions, as introduced below in Sec. 1.1, 
was first studied by J. Ambjorn and J. Jurkiewicz back in 1992 [1], They found that the model 
possesses two phases and the transition between them initially seemed to be of 2 nd order, which is 
necessary for a continuum limit to be defined. In 1996 then, P. Bialas, Z. Burda, A. Krzywicki and 
B. Petersson reported for the first time the finding of some 1 st order behavior in this phase transition 
for systems consisting of IV 4 = 32k 4-simplices [5], Shortly afterwards B. V. de Bakker verified this 
finding and extended the study to larger systems with A 4 = 64k. However, we were not completely 
convinced by the numerical methods used in the latter work. In particular, there were three things 
which disturbed us: 

1. Measurements were taken after a fixed number of accepted (instead of attempted) moves, 
which introduces a systematic error. 

2. The use of an artificial harmonic potential to control volume fluctuations also introduces a 
systematic error. 

3. Autocorrelation and thermalization times could easily have been underestimated. 

Therefore we wanted to cross-check these old results with our own, hopefully correct code which 
satisfies detailed balance, uses a potential well instead of a harmonic potential to control volume 
fluctuations, and makes use of parallel tempering to cope with critical slowing down. 

The paper is organized as follows: 

In the remainder of this section we give a brief overview of the EDT model and its phase diagram. In 
section 2 we describe our simulation methods while in section 3 we present our results: after having 
verified in part 3.1 that the phase transition is 1 st order, we address in part 3.2 the question whether 
we can observe a coexistence of the two phases and present therefore a local criterion to determine 
whether a piece of triangulation is in a crumpled or elongated state and try to identify the nature of 
the metastability in the Markov chain that causes the 1 st order transition. Finally in the appendix, we 
propose a modification of the path-integral measure, based on a counting of the number of possible 
moves in each triangulation, which could weaken the 1 st order nature of the phase transition. 


1.1 The EDT Model 


In 4-dimensional Euclidean Dynamical Triangulation (EDT) [ i ] the formal path integral for Eu¬ 
clidean (local SO (4) instead of SO (3,1) symmetry) gravity, 


Z = 



( 1 . 1 ) 


with the Einstein-Hilbert action Seh'- 

Seh = ~ iq^q J d4 x Vd {R - 2 A), 


( 1 . 2 ) 


is regularized by approximating the configuration space (space of all diffeomorphism inequivalent 
4-metrics) with the space of simplicial piecewise linear (PL) manifolds consisting of equilateral 
4-simplices with fixed edge length a (such manifolds are also called abstract triangulations). Un¬ 
der such a discretization, (1.2) turns into the Einstein-Regge action Ser which for equilateral 4- 
simplices and a space-time of topology S 4 takes the simple form 


Ser = —K 2 N 2 + K 4 N 4 . 


(1.3) 


Ni labels the number of /'-simpI ices in the PL manifold and re 2 , mi are related to the bare gravitational 
and cosmological couplings G, A by 


k 2 = 


V2_ 
4 G 


K4 


10 arccos(l/4) V 2 + AV 4 
87 tG 


(1.4) 
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with V n = a n being the volume of a n-simplex. The partition function (1.1) can now be 

written as 

Z(k 2 ,« 4 ) =Y,7j- = ^2 Z{k 2 ,N±) e~ KiNi , (1.5) 

T T N 4 

where after the first equality sign the sum runs over all abstract triangulations T of S 4 and Cp is a 
symmetry or degeneracy factor (to avoid over-counting) which is assumed to be ~ 1 for sufficiently 
large systems. After the second equality sign, the canonical partition function 

Z(k 2 ,N 4 ) = Y. ( 1 . 6 ) 

{T:JV 4 (T)=Af 4 } L ' T 


was used. 

The partition function (1.5) is suitable for use in a Markov chain Monte Carlo simulation with 
Metropolis updates consisting of the so-called Pachner moves (see Sec. 2.1). 


1.2 Phase Diagram 

The grand canonical partition function (1.5) is finite only if K 4 > nf (n 2 ). We therefore have a 
critical line for convergence in the (k 2 , K 4 )-plane, given by nff (n 2 ). To obtain the thermodynamic 
limit (JV 4 —» 00 ) we have to ensure that K 4 (k 2 ). For quasi-canonical simulations around 

some fixed volume A r 4 1 , we can use (1.6) to define a pseudo-critical K p f r by 


«r(tc 2 , n a ) 


5In (Z(k 2 , N 4 )) 
dN 4 


N4—N4 


(1.7) 


which corresponds to the value of k 4 for which the ^-distribution is flat around N 4 . 



Figure 1: Phase diagram for 4D EDT. The figure shows K pcr (K 2 , N4) as a function of K 2 for different N4 
together with the corresponding pseudo-critical points («f cr (JV4), K pcr (N 4) , IV4)). The dotted red line 
separates the crumpled from the elongated phase; in the limit N4 — > 00 this line ends at the critical point'. 
(k 2 t , K| r ). To improve readability, the y-axis shows (k 4 — 2.33 k 2 ) instead of K 4 itself. 

For constant JV 4 , we can define (see Fig. 1) a line K pcr (k 2 , N 4 ) as a function of k 2 , along which 
two phases are separated by a pseudo-critical point at n 2 = n^ cr (N 4 ). For k 2 < (N 4 ) we 

are in the crumpled phase where a typical configuration is highly collapsed in the sense that the 

1 There is no set of ergodic moves known for fixed volumes, it is therefore necessary to let the volume fluctuate around the 
desired value N4. 
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distance between any two 4-simplices is very short, leading to a large (infinite) Hausdorff dimen¬ 
sion. For K2 > ^2° (-V4) we are i n the elongated phase with Hausdorff dimension ~ 2 , where a 
typical configuration consists of a so-called baby-universe tree: the total volume is subdivided into 
smaller parts, the baby-universes * 2 , which are pairwise connected by only a small minimal neck 3 . 
This structure is hierarchical in a treelike manner: consider the largest baby-universe as “mother” 
with outgrowing smaller “babies” which in turn give birth to their own “babies”, and so on (see 
right-hand side of Fig. 2 ). 

The true critical point in the thermodynamic limit is obtained as 

(«?, «D = lim («T (iV 4 ), «T («r m , N 4 )) . (1.8) 

N 4 — >00 



Figure 2: Representative configurations in the crumpled (left, K 2 = 1.26) and elongated (right, K 2 = 1.30) 
phase at system size N 4 « 64k: in the crumpled phase, the triangulation consists of one large, highly connected 
bunch with outgrowths which are at least an order of magnitude smaller. In the elongated phase on the other 
hand, although a largest component still exists and may be called “mother universe”, it is much smaller than in 
the crumpled phase and some of its outgrowths (the “babies”) are of comparable size. 


2 Simulation Methods 

2.1 Pachner Moves 

In d dimensions, there exist (d + 1) Pachner moves. They form an ergodic set of local updates 
[15] in the space of abstract triangulations of fixed topology without boundary. A n-move (n £ 
{0 ,..., d}) consists of the following procedure: pick a n-simplex which is contained in (d + 1 — n) 
d-simplices. The complex consisting of these (d + 1 — n) d-simplices has the same boundary as a 
corresponding dual complex spanned by (n + 1) d-simplices that share a common (d — n)-simplex. 
We can therefore just replace the complex around the selected n-simplex with its dual (see Fig. 

3 for an illustration in four dimensions). The only additional constraint is the so called manifold 
constraint , that is: the (d — n)-simplex shared by the (n + 1 ) newly created d-simplices of the dual 
complex must not already exist in the triangulation as this could result in topology changes. From 
now on we will consider only the 4-dimensional case. 

-We consider a baby-universe as a collection of 4 -simplices which could all be pairwise connected by a path on the dual 

lattice that does not pass through any 3 -simplex that belongs to a minimal neck. 

3 In four dimensions, a minimal neck consists of five 3 -simplices forming a closed hyper-surface that looks like the bound¬ 
ary of a 4 -simplex, but without a corresponding 4 -simplex being present in the triangulation. Intuitively, a minimal neck is 
something like the bottleneck of a sand glass (but an extremely narrow bottle neck). 
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Figure 3: Illustration of a Pachner 3-move in four dimensions: on the left-hand side we have two 4-simplices 
which share a common 3-simplex (red, projected onto two dimensions). The 3-move replaces this complex with 
the one shown on the right-hand side of the figure where four 4-simplices share a common 1-simplex (blue). 
The 1-move is just the inverse of the 3-move. 


2.2 Detailed Balance 

Calling T k the current triangulation in our Markov chain, we obtain T k+ \ as follows: 

1. randomly choose a move type n £ {0,..., 4}, 

2. randomly choose one of the A4 4-simplices of 7), and call it D, 

3. randomly choose one of the ( n ^ 1 ) n-simplices contained in D and call it S , 

4. perform a Metropolis test with acceptance probability p n (T k , S ): 

• accept: T k+ \ is obtained from T k by applying the n-move at S, 

• reject: T k+1 = T k . 

The acceptance probability at step 4 is given by [ 10] 


Pn ( T,S) = 


[ p n (7V 4 (T)) if n-move possible at S G T 


else 


( 2 . 1 ) 


where p n (N4) = min jl, jv 4 +AJV 4 (ri.) e KaAArK 4 Aiv 4 (n) j j s so . ca jj ec i reduced transition 

probability, AN, (n) labels the change of N, under a n-move, and a n-move is considered as possible 
at S' if S' is contained in (5 — n) 4-simplices and the application of the move does not violate the 
manifold constraint mentioned above in Sec. 2.1. 

Equation (2.1) can be derived from the detailed balance condition: assume first we have two valid 
4-dimensional triangulations T , T' , where T' can be obtained from T by applying a n-move at a 
specific n-simplex S of T. The detailed balance equation then reads 


p(T) P (TAT') = p (T') P (V A A n T) , 


( 2 . 2 ) 


where p (T) = e K2 N2 ^ 7 "> K4 N A T ) anc j p (T AT') is the transition probability. The latter can be 
written as 


P (TAT') = 


5 — n 


5 N± (T) ( n 5 +1 ) 


Pn (N 4 (T)) , 


(2.3) 


where p n (N 4 ) is again the reduced transition probability and the factor in front of it is the probability 
for selecting (with the update scheme mentioned above) the n-simplex S through which T and T' 
can be related by applying a n-move: 1/5 is the probability for choosing the correct move type n, 

r ‘r/' the one for selecting a 4-simplex D which contains S and , J . is the probability for selecting S 

4 U+i) 

Qp rhe f 1 ^_pimnli/'oc /-.f D ^ n fU*. n rr,_cimnlAvt 1 


n-simplices of D. Note that 


is the local 4-volume of a n-simplex 4 that allows 


V„+1, “*■*“*'*'» —^ CT 

for a n-move (in units of 14) and p n (N 4 ) can therefore be interpreted as the transition probability 


4 The 4 -volume containing all points which are closer to the n-simplex under consideration than to any other n-simplex. 
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for a n-move per unit volume, as the term reduced suggests. 

For the inverse transition probability on the right hand side of (2.2) we have analogously: 


P(T' 4 A"T) = 


n + 1 


5 ^4 (TO ( 5 !J 


PA—n (N 4 (T')) = 


5 — n 


5 ^4 (TO (4J 


P4-n (iV 4 (T')) • (2.4) 


Equation (2.2) reduces therefore to 

P ( T ) mr) Pn (Na (T)) = P (TO jv^Ty PA-n (JV 4 (TO) , (2.5) 

which, by noting that JV4 (TO = AO (T) + A AO (n), is satisfied by setting 

= mi "{ 1 ' JV 4 + AiV 4 (n) e " A1V,( ’ ,l ~* ,A ' , ‘ < " ) }- (2 ' 6) 

This gives the upper part of (2.1). For the lower part, when T' does not exist, we have p (T') = 0 
and in order to satisfy ( 2 . 2 ), we have to set 

P(TAT') = - 5 ~ n / . , p n (T, S) = 0. (2.7) 

5 AO (T) ( n 5 +1 f nK 


We now have a prescription for how to produce a Markov chain containing the configurations 
required to evaluate (1.5). As neighboring elements in our Markov chain are highly correlated, it is 
appropriate to take measurements only on every fc th element in the chain, where k must be a constant 
in order to preserve the probability distribution. As already mentioned at the beginning, this was 
not always respected in previous work, as e.g. in [9] the measurements were separated by a fixed 
number of accepted moves, which turns the true separation k between measurements into a random 
variable whose value depends on what kind of configurations is currently sampled by the Markov 
chain. 


2.3 Controlling the Volume 

As K pcr (« 2 , N 4 ) is monotonically growing with N4, it is practically impossible to run fully grand 
canonical MC simulations for the EDT model. But canonical simulations are anyway better suited 
to investigate finite size scaling. Unfortunately, as already mentioned before, there is no set of 
ergodic moves known for the space of triangulations of fixed N4 and it is therefore not possible to 
run canonical simulations either. The best we can do is to run quasi canonical simulations based on 
(1.5) but with N4 constrained to fluctuate around some desired N 4 . In previous work [1,5,9], this 
was often achieved by adding a harmonic potential, 

U(N 4 ,N a ,S) = 6 - (N 4 - N4) 2 , (2.8) 

to the action (1.3). This of course introduces a systematic error for all moves which change N4. We 
therefore decided to rather use a infinite potential well of some reasonable width w ~ 2 a (N 2 ) / 2.5, 

where 2.5 = max | ^^{n) } anc ^ a (-V 2 ) is the square root of the ^-variance. 

As with such a potential well we cannot use the saddle point expansion method from [10] to tune 
K 4 to its pseudo-critical value n P 4 r (k, 2 , N 4 ), we instead made use of a method mentioned in [4]: as 
the A^-histogram has to be flat around N4 if K4 = ( K 2 , N 4 ), we have that 

pT° (Na) pT (N 4 ) = P 9 0 eo (N 4 + a N 4 (4)) P p 0 cr (N 4 + A N 4 (4)), (2.9) 

where pPp r (N4) is the reduced transition probability (2.6) with K4 = K pcr (k 2 , N 4) and j?® eo (N4) 
is the average geometric probability for a n-move, i.e. the fraction of the total volume that allows 
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for a change (through one of the Pachner moves), to which a n-move could be applied 5 (see Fig. 4): 


P 9 n°(T) = 


(r) 


( 2 . 10 ) 


E Nm (^) f™ 9al (T) 

n =0 \m+l) 


m —0 


where fn 9al (T) is the fraction of n-simplices in the triangulation T where a n-move is possible. 
One can then solve for n pcr (« 2 ,-/V 4 ) which leads to 


K 


r (« 2 ,jv 4 ) = 

1 


A N 4 (4) 


In 


pT° (n 4 


vl eo (N 4 + A N 4 (4)) 


— In ( 1 + 


A N 4 (4) 
N 4 


A N 2 (4) 
A N 4 (4) 


« 2 . ( 2 . 11 ) 


Note that the 4-move is always possible and one therefore effectively only has to measure the average 
fraction of vertices which allow for a 0-move in order to determine K pir (kv. N 4 ). Furthermore, as 
the average geometric probabilities vary slowly with N 4 , we can, as long as the width of the potential 

instead of =— 


well for N 4 is much smaller than N 4 , just use 
pseudo-critical value. 


(pc 9 E 


p a 0 co (N i +AN i ( 4 )) 


to set K 4 to its 


2.4 Autocorrelation Time 

The autocorrelation time r is a measure for the typical distance between two uncorrelated elements 
in a Markov chain. It can be thought of as the time it takes for a change to propagate through the 
typical volume of the system over which the degrees of freedom are correlated. One therefore writes 
r oc £ z , where £ is the correlation length and the dynamical critical exponent z is expected to be 
z ~ 2 for a local update scheme; i.e. updated information propagates in a diffusion-like manner. 

For a 2 nd order transition where the correlation length diverges when approaching the critical 
point, £ is truncated by the linear system size L and we have r ex L z = V Z/ ' d " with dn being 
the Hausdorff dimension of the system. In this case the autocorrelation time obviously diverges as a 
power of the system size. 


For a 1 st order transition £ remains finite at the transition point for all system sizes. Nevertheless 
the autocorrelation time can diverge even more dramatically as transitions between the two phases 
that coexist at the transition point become exponentially suppressed with increasing system size. The 
autocorrelation function should then consist of two parts: a relatively steep first one corresponding 
to the decay of autocorrelations within a single phase, as well as a second, much less steep part 
which reflects the fact that the system remains in one and the same phase for a rather long time. 
Unfortunately it is almost impossible to verify this as it would need ridiculously long simulations to 
obtain the required accuracy on the auto-correlation function. 


In both cases, for 1 st and 2 nd order phase transitions, parallel tempering can be used to reduce 
autocorrelations [11,12]. The idea is to run K simulations for different values of the couplings in 
parallel where the couplings are chosen such that they lie on a line in coupling space which connects 
a region with slow relaxation with another where relaxation is fast. One now periodically attempts to 
swap configurations between neighboring sets (called replicas) with an acceptance probability given 


'’Alternatively, one could define p® e ° (T) as the fraction of the overall total volume of a triangulation T, to which a 71- 
move could be applied, which would lead to a different normalization: p® e ° (T) = r> N n ( T ) f^ 9al (T), such 

that p® eo (T) would coincide with the probability to select a good location for a n-move within the update scheme described 
above in Sec. 2.2. 
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k 2 


Figure 4: The average geometric probabilities p® eo (IV4), n £ {0,..., 4}, are plotted as functions of «2 for 
different system sizes N 4 = 2k (blue), 4k (red), 8 k (dark gray), 16k (dark yellow), 32k (dark blue), 48k (dark 
red), 64k (dark green). Note that, as the 4-move is always possible, the change in pj[ eo ( N 4 ) is just due to the 
normalization chosen in ( 2 . 10 ). 

The last plot at bottom right shows the pseudo-critical K% cr (k 2 ,N 4 J corresponding to p q E ° (-/V 4 J and 
pf° (N 4) as given by (2.11). It is a close up version of Fig. 1, showing in more detail the behavior of 
k v cr (ft 2 , N 4 ) close to the pseudo critical point (to improve readability, the y-axis shows again (t ?4 — 2 . 33 /t 2 ) 
instead of K 4 itself). The small red dots indicate the pseudo-critical points. 
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by 


p SWap [{(«*,iV 4 ) , («[,JV')} ->■ {(^X) , («<,JVi)}] 

= min{l J e^’ Jv ^- s ( K - JV 0+ s «- JV 0- s «’ JV 0} I (2.12) 

where (k, : . AT), («', A 7 ') are the sets of couplings and configuration variables for the two neighbor¬ 
ing replicas and S (re*, AT) is the action of a configuration with variables N, at couplings k. l . The 
advantage of this procedure is twofold: first, we no longer have just one Markov chain per simulation 
point in coupling space that has to build up the whole corresponding statistical ensemble but now 
all K chains alternately contribute to all of the statistical ensembles at different couplings. Second, 
if every Markov chain frequently reaches regions in coupling space where relaxation is fast before 
passing again through the critical region, then the configurations which the chain contributes to the 
statistical ensembles in the critical region, are on average much less correlated than corresponding 
configurations of a Markov chain that remains all the time in the critical region. For more details 
see [11,12] where it is also explained how this procedure can be optimized. Especially in [12] the 
application of parallel tempering to 1 st order transitions is discussed. 

In our implementation we chose, for a fixed average volume AT, 24 or 48 equally spaced (w.r.t. 
the k 2 direction) couplings along the pseudo-critical line K pcr (n 2 - N 4 ), such that they join a region 
in the crumpled phase where relaxation is fast, with one in the elongated phase where relaxation 
is also fast (compared to the critical region), and thereby pass through the critical region around 
K% cr (AT). After some runtime, the optimization procedure of [ 1 1 , 12 ] is applied which gives us a 
new set of couplings for which the configuration exchange between replicas is more frequent. 

In Fig. 5, two Monte Carlo time histories for the observable AT (number of triangles) are shown 
for comparison, both for a system of size N 4 = 32k at the pseudo-critical point k 2 = K pcr (32k) « 
1.258. The left one stems from an ordinary simulation while the right one was obtained using parallel 
tempering. 
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Figure 5: The figures show Monte Carlo time histories for the observable N 2 (number of triangles) for systems 
of size AT = 32fc with K 2 set to the pseudo-critical value K pcr (32k) ss 1.258. The left-hand figure corresponds 
to the result of an ordinary simulation whereas the right-hand figure was obtained using parallel tempering with 
24 replicas. 


2.5 Data Analysis 

Due to the use of a potential well instead of a harmonic potential to control the system volume and 
due to the tuning of K 4 to its pseudo-critical value, we have significant volume fluctuations in the 
data which also affect for example the AT distribution. To take this into account, we project the data 
in the (AT, AT )-plane along the "correlation direction" before evaluating any observables (see Fig. 
6 ), i.e. instead of AT we use 

N 2 = N 2 -f(N 4 ) (2.13) 
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to evaluate observables depending on N 2 , where 



(N 4 - (iV 4 )). 


(2.14) 


We checked that this leads to the same results as when evaluating the observables only on data sub¬ 
sets corresponding to single, fixed N 4 values. 
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Figure 6 : The figure shows the N 2 -N 4 distribution at k 2 = 2.8014 for a system of average size IV 4 = 64k. 


As long as the fluctuations in N 4 are forced to be much smaller than the average system size IV 4 itself, we can 
project the data along the A^-AT-correlation direction (indicated by a red line) and evaluate observables as if 
we had a true canonical simulation at system size A 4 and fluctuating triangle number N 2 as given by (2.13) 
and (2.14). 

After that, we use multi-histogram reweighting [13] with respect to n 2 . The parallel tempering 
optimization procedure mentioned above also leads to a good distribution of simulation points for the 
reweighting. The errors are determined by the Jack-Knife method with 20 sets. In multi-histogram 
reweighting, these sets consist of the simultaneous data of all the simulations at different n 2 values 
(i.e. to form the Jack-Knife sets we consider as a measurement all the measurements at different 
k 2 values which correspond to the same Monte Carlo time), therefore cross correlations should 
automatically be taken into account. 

3 Results 

3.1 Order of Phase Transition 

3.1.1 N 2 Distribution 

As stated in [9], the pseudo-critical ^-distribution starts to be visibly double-Gaussian for systems 
consisting of more than about 32000 4-simplices. For systems containing 32k, 48k and 64k 4- 
simplices, these distributions are shown in Figure 7, where k 2 was tuned to produce peaks of equal 
height (left) or equal area (right). Either way it can be seen that the double peak structure becomes 
more pronounced with increasing system size and that there is so far no sign that the peaks will 
merge again in the thermodynamic limit. This behavior is characteristic of a 1 st order transition. 
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Figure 7: Normalized ./V 2 -distribution for systems of size N 4 = 32k (blue), 48k (red) and 64k (green): the 
solid lines are double-Gaussian fits to the data. To the left, the values of re 2 were chosen such that the two 
Gaussian parts of the distribution function have the same height, whereas on the right-hand side the «2 values 
are such that the two Gaussians have the same area, i.e. the two states are equally probable. It can be seen that 
the double peak structure becomes more pronounced with increasing system size and there is no sign that the 
peaks will merge again in the thermodynamic limit. This is characteristic of a 1 st order transition. 


The reason for the double peak structure in the N 2 -distribution is, that at a 1 st order transition 
point, the two phases can coexist (see Fig. 8 ) while N 2 takes on different average values in each of 
these phases. 


-¥ 



Figure 8 : Example of a crumpled (left) and an elongated (right) configuration with IV 4 m 64k, both recorded at 
re 2 = 1.28, which is approximately the pseudo-critical point. 


3.1.2 Scaling of f? 4 


A more quantitative method to determine the order of a phase transition is to study finite-size scal¬ 


ing of the 4 th order Binder cumulant (Kurtosis) of the N 2 distribution, /i 4 [A)>] = 
According to [14], for large enough N 4 this quantity should scale like 


{(N 2 -(N 2 )) i ) 


(( N 2 -{N 2 )) 2 ) 2 ’ 


BT [iV 2 ] (JV 4 ) « Bf [N 2 ] + d jV 4 -", 


(3.1) 
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where _B| r [iV 2 ] is the critical, infinite volume value of the Binder cumulant. For a 2 nd order transition 
one should get 1 < B"' [A/ 2 ] < 3 and u = 1/ dnu, where v is the critical exponent of the correlation 
length £at 2 ~ {nf — k 2 | v and dn the Hausdorff dimension, whereas for a 1 st order transition we 
should obtain [A/ 2 ] = 1 and u> = 1. 

We tried to fit our data for B% cr [A T 2 ] (A/ 4 ) assuming 1 st and 2 nd order scaling ansaetze (see figures 
10 and 9). For the 2 nd order ansatz, we looked at B^ cr as a function of the average linear system size 
Lpcr i ns tead of volume N 4 , as a diverging correlation length will be truncated to this L pcr , which is 
defined as 

L pcr (JV 4 ) =^r«(r; JV 4 , W), (3.2) 


with n (r; A/ 4 , n 2 ) being the average volume profile of a triangulation of size N 4 at a given value of 
k 2 , i.e.: 


(r; A/ 4 , k 2 ) = ( 1^2 H d A(s!,S 2 )-r) 


Nl 


(3.3) 


where si, S 2 run over all N 4 4-simplices and d& (si, s 2 ) is the geodesic distance between the sim- 
plices si and s 2 in the triangulation A and (...) refers to the average over triangulations at «: 2 . As 
is typical for a weak 1 st order transition, the 2 nd order fit seems to work fine, but the obtained values 
B^ [N 2 ] = —4.2 ± 5.0 and v = = 2.01 ± 0.74 do not make much sense. Instead fixing 

u = 1 for the 1 st order ansatz, and using only the data points from the largest two simulated systems 
(A /4 = 64k and 48k), we obtain Bf [A/ 2 ] = 1.00 ± 0.08 which is the expected value for a 1 st order 
transition. 



Li” 


Figure 9: Binder cumulant B^ cr [JV 2 ] as a function 
of average linear system size L pcr (N 4 ) (3.2), assum¬ 
ing a 2 nd order transition, together with a fit of the 
form (3.1) where N£ = ^P cr 'j 1 / L ' _ \y e a j so in¬ 
cluded higher order corrections. It can be seen that 
the fit seems to work fine, but the obtained values 
B% r [N 2 ] = -4.2 ± 5.0 and v = 2.01 ± 0.74 
do not make much sense. 


64k48k 38k 16k 8k 
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Figure 10: Binder cumulant B pcr [A/ 2 ] as a function 
of I/N 4 . The red lines correspond to fits of the form 
(3.1)withtJ = 1 (1 st order transition) to the data of the 
largest and second largest pair of systems. The value 
B1 r [A/ 2 ] = 1.00 ± 0.08 obtained from the largest 
pair is compatible with the expected value of 1 for a 
1 st order transition. The fit for the second largest pair 
yields a too large value: [A/ 2 ] 1.50 ± 0.07. This 

large change indicates that we are still far away from 
the thermodynamic limit. 


A fit to the data of the next smaller pair of systems, i.e. those consisting of 48k and 32k 4- 
simplices, leads to _Bf r [A/ 2 ] = 1.50 ± 0.07. This large change in the value of the extrapolated 
/f]' r [A/ 2 ] indicates that we are still far away from the thermodynamic limit. 
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3.2 Coexistence of Phases 


As the phase transition is 1 st order we expect for finite systems a coexistence of the elongated and 
crumpled phases in some neighborhood of the pseudo-critical point. Instead of speaking of crumpled 
and elongated configurations as in Fig. 8, we should rather speak of configurations in which the 
crumpled or the elongated part dominates. For example, the left-hand graph in Fig. 8 is dominated 
by the large bubble 6 in the middle that is in the crumpled state, but attached to that large bubble are 
baby-universe trees which are in the elongated state. Similarly, the right-hand graph in Fig. 8, whose 
strong branching indicates that the configuration is mainly in the elongated state, also contains some 
larger bubbles that presumably are in the crumpled state. 

Although criteria like "strong branching" and "large bubbles" seem to work well to decide if a piece 
of triangulation corresponds to the elongated or crumpled phase, there is some ambiguity in what 
"strong branching" should mean or what bubble sizes should be considered as small. As a first 
attempt, we could define the elongated phase as consisting of bubbles of size 6 and everything else 
as belonging to the crumpled phase. By the "size" of a bubble, we mean from now on the number 
of the bubble’s 4-simplices plus the number of its minimal necks 7 (i.e. the volume the bubble would 
have after replacing all minimal necks by ordinary 4-simplces). In contrast: the "volume" of a bubble 
still refers to just the number of 4-simplices of that bubble. Figure 11 shows, as a function of k 2 , 
the average fractional volume of a triangulation contained in size 6 bubbles. The complementary 
fractional volume would therefore be the one contained in bubbles of size larger than 6. It can be 
seen that at small values of fc 2 , deep in the crumpled phase, the volume is dominated by contributions 
from large bubbles (in fact one very large bubble containing almost all the 4-simplices) whereas at 
large values of k 2 , deep in the elongated phase, the dominant contribution to the total volume comes 
from the size 6 bubbles. The swap in the dominance occurs not exactly at the pseudo-critical point 
but at a somewhat larger value of fc 2 = fc? which, however, seems to coincide with the expected 
infinite volume critical value of the coupling nf « 1.33. 

As the curves in Fig. 11 show no volume dependency in the elongated phase and may get close to 




Figure 11: The left-hand figure shows the ensemble averaged fractional volume contained in bubbles consisting 
of five or less 4-simplices (i.e. bubbles of size 6 ) as a function of fc 2 at system size N 4 = 32k. With increasing 
ft 2 , it can be seen that at some point K 2 = ft? (IV 4 ), the curve goes above 0.5 and therefore the two classes 
of bubbles, those consisting of more than five and those consisting of five or less 4-simplices, change their 
roles as dominant and non-dominant parts of the system. The right-hand figure shows the same quantity but for 
three different system sizes, N 4 = 32k (dark blue), 48k (dark red), 64k (dark green), in a neighbourhood of 
the corresponding pseudo-critical points: ft 2 cr ({32k, 48k, 64k}) = {1.258,1.271,1.280}. Due to the finite 
volumes, ft 2 cr (IV 4 ) is smaller than ft? (IV 4 ) which seems to be volume independent and is very close to the 
expected infinite volume critical value of the coupling ft 2 r ss 1.33. 


unity only in the limit k 2 —¥ 00 , it should be clear that the characterisation of this phase as consisting 
of just size 6 bubbles is not adequate. The reason is, that if a system would consist of size 6 bubbles 
only, Pachner 3-moves could only be applied to 3-simplices which are part of minimal necks, such 
that these moves would necessarily destroy those necks and thereby produce bubbles of size larger 

6 Bubble is just another word for baby-universe as defined above in footnote 2. 

7 A bubble of size 6 (which is the smallest possible size) can consist of five 4-simplices plus a minimal neck or four 
4-simplices plus two minimal necks, and so on. 
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than 6. As can be seen in Fig. 12, which shows as a function of re 2 the quantity A necks (n), i.e. the 
average change of the number of necks under a n-move s , indeed, more and more 3-moves change 
the number of necks as k 2 increases. But nevertheless, the size 6 bubbles seem to be the dominant 
building blocks of the elongated phase. 



K 2 


Figure 12: The figure shows for a system of size N 4 = 32k, as a function of k 2 , the average change of the 
number of necks caused by the next n-move. The 0- and 4-moves always remove or add a "volume 5" bubble 
and a corresponding neck, while the 2-move never changes the number of necks. For the 1- and 3-moves 
A necks (n) changes as function of k 2 as the fraction of 3-simplices, which allow for a 3-move and are also 
part of a minimal neck changes. For k 2 > 1.6, it seems that the average 3-simplex which allows for a move, is 
rather part of more than one neck than of no neck, which is why A necks (3) drops below -1. The behaviour 
of A necks (1) follows from the fact that the 1-moves is the inverse of the 3-move. Thus for k 2 > 1.6, the 
triangle that will be created by applying a 1-move to one of the 1-simplices that allow for such a move, will 
be rather part of more than one neck than of no neck. At the pseudo-critical point = 1.258, we have that 
A necks ({1, 3}) « {0.8, -0.8}. 

In the pseudo-critical region, in order to change from a rather crumpled to a rather elongated 
state, the system has to go through the process of producing and growing new baby-universe branches 
on top of the large "mother universe", until almost the whole volume fits into them while a (distin¬ 
guishable) "mother universe" disappears. This is illustrated in Fig. 13 where, as a function of k 2 , it 
is shown how the total volume of a system with N 4 = 32k is distributed, on average, over bubbles 
of different sizes. 

Regarding this process of creating new baby-universe branches, which necessarily starts with the 
creation of a new size 6 bubble, it is interesting to note that applying a 4-move to a 4-simplex con¬ 
tained in a bubble that consists of only five 4-simplices results in a triangulation which is slightly 
more restrictive with respect to further application of 3- and 0-moves, as compared to the case where 
a 4-move is applied outside of such a "volume 5"-bubble. This is explained in more detail in Fig. 14 
for the two-dimensional case, and in figure 15, we show that the effect of this mechanism is indeed 
observable: the latter figure shows, for different system sizes, the average numbers of possible Pach- 
ner n-moves at the pseudo-critical point as a function of N 2 /N 4 . Comparing the graphs in Fig. 15 
with Fig. 7 in order to identify which A^/A^-interval corresponds to which phase, we see that the 
numbers of possible moves undergo an abrupt change precisely in the region where the valleys of the 
corresponding graphs in Fig. 7 are located. These abrupt changes are just as one would expect from 
the phenomenon described in Fig. 14 which occurs a soon as a 3-move is applied to a 3-simplex that 
is part of a minimal neck between a "volume 4"- and a "volume 5"-bubble (which leads to a "volume 
ll"-bubble i.e. a "size 12"-bubble with just one neck): the number of possible 3- and 0-moves is 
lower and the number of possible 2-moves higher than in a configuration of the same size but without 
this particular "volume 11 "-bubble. 

At the beginning of this section, we mentioned also using the branching factor, at least on an 

s The quantity A necks ( n ) is computed by counting for each location where a n-move could be applied, the number of 
necks that would be created or destroyed by such a move, and taking the average. 
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Figure 13: The figures both show (on different scales) as a function of K 2 how the total volume of a N 4 = 32k 
system is distributed, on average, over bubbles of different sizes, where the size is given by the number of 4- 
simplices plus the number of minimal necks: for k ,2 « k ! ^ t (-/V 4 ), almost the whole volume is concentrated 
in just one large bubble, but with increasing K 2 the volume distributes over more and more (and larger) small 
bubbles until the (distinguishably) largest bubble disappears for K 2 ~ «;f cr ((V 4 ). At this point, the largest 
bubble contains only about 20 % of the total volume. 


intuitive level, as a criterion to decide if a piece of triangulation is in the elongated or crumpled 
phase. The branching factor itself turns out not to be a good criterion to distinguish between the 
two phases as its average value drops again with increasing K 2 for K 2 > n^ cr (see Fig. 16). On the 
other hand, the related average neck density or "branching factor per size" is monotonic and seems 
to yield a meaningful criterion to distinguish between the two phases (see Fig. 17). In figure 18 we 
show how the total volume distributes over bubbles with different neck densities but it remains to be 
understood what the exact value of the critical density, which seems to be around 1/9, should be and 
where it comes from. 
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Figure 14: Illustration in two dimensions: in the left Figure, by applying a 2-move to the general triangle 
(abc) (which is not part of a "volume 3"-baby-universe), we have replaced the triangle by a baby-universe "A" 
(horizontally shaded) consisting of the three triangles (abg), (beg) and (cag). We could continue by applying 
1-moves to all the 1-simplices of the neck (abc) which would replace the link (ac) by (gf), (cb) by (gd) and (ba) 
by (ge). In the right hand figure, we have instead produced another "volume 3 "-baby-universe "B" (vertically 
shaded) with neck (age) inside the already existing baby-universe "A". In this case, we can only apply a 1-move 
either to (ac) and (ag) or to (ac) and (gc), as applying a 1-move to both (ag) and (gc) would result in a double 
link (hb). 

It is also clear that the insertion of the vertex "h" into one of the triangles of "A" makes it impossible to remove 
the vertex "g" by a 0-move. If "h" had instead been inserted into e.g. the triangle (acf), not just "h", but also "g" 
could still be removed. 

Similarly in the four dimensional case, a piece of triangulation produced by applying a 4-move inside a "volume 
5"-baby-universe "A" to produce a new "volume 5"-baby "B” of "A", has the following effects: 

- the number of possible 0-moves is not increased, as the 4-move that created "B" has also destroyed an 
already existing location where a 0-move could have been applied: the vertex that was in the center of 
"A". This does not happen, if a 4-move is applied to a 4-simplex that is not part of a "volume-5"-bubble, 

- all the 3-simplices that are only part of the neck between "A" and "B" (but not of the neck between "A" 
and the rest of the triangulation), would allow for a 3-move, but by performing one of these 3-moves, 
the remaining ones become impossible. If "B" were not the baby of a "size 5,4 or 3"-bubble, 3-moves 
could in general be applied to all of its neck’s 3-simplices, 

- and finally, in four dimensions, there is also an effect on the 2-move: the application of a 3-move to one 
of the 3-simplices that are part of only the neck between "A" and "B", removes another location where a 
0-move had been possible (as it destroys "B"), but at the same time, generates three new locations where 
a 2-move could be applied. 
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Figure 15: From top to bottom: average numbers of possible 3,2,1 and 0-moves (normalized by the system 
size) vs. N 2 /N 4 for systems of size (from left to right) N 4 = 32k, 64k and 128k at the pseudo-critical point. 
Comparison with Fig. 7 shows, that the strong changes (.e.g. for N 2 /N 4 = [2.398, 2.402] in the N 4 = 64k 
case) happen at the location of the valley of the corresponding graph in Fig. 7. 
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Figure 16: The figure shows, as a function of K 2 , the 
average branching factor, i.e. the average number of 
necks, of bubbles which are neither the largest bubble 
in the system nor volume 5 bubbles, which are just the 
terminating leaves of a baby-universe branch. As the 
branching factor decreases again after the phase tran¬ 
sition at K 2 « 1.258, a large branching factor alone is 
not a good indicator for a bubble to be in the elongated 
phase. 


Av. number of Necks of 'Motheruniverse' 



Figure 17: The figure shows, as a function of K2, the 
average density of necks of the largest bubble in the 
system, i.e. "number of necks of bubble"/"size of bub¬ 
ble" (where the size is again given be the sum of the 
numbers of necks and 4-simplices). As the largest 
bubble can be assumed to correspond to the crumpled 
phase for K 2 < n%f r (IV 4 ) « 1.258, but for K 2 » 
(A 7 i) is just the slightly largest of many almost 
equally sized bubbles, which are all part of the elon¬ 
gated phase, this shows, that it is the neck-density of a 
bubble rather than its total neck number which distin¬ 
guishes between the two phases. Note also that for size 
6 bubbles, the neck density is always > 1/6 ~ 0.167 
which according to this figure is clearly elongated, as 
it should be. 



Figure 18: The figure shows for a system of size IV 4 = 32k, how the total volume is distributed over bubbles 
with different neck densities, i.e. different ratios "number of necks of the bubble'Vsize of the bubble" and how 
this distribution changes with K 2 . 
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3.3 Balls in Boxes Model 


It is well known that the balls in boxes model [ 6 , 7] describes nicely the qualitative features observed 
in EDT simulations. In the canonical formulation the model describes the statistical ensemble of a 
fixed total number N of balls distributed in a varying number M of boxes: 

OO 

Z(N,k) = Y e KM Y P(li)-"P( < lM)SN,q 1 +...+ qM , (3.4) 

M=1 q\ 

where p (q) is the probability for a single box to contain q balls and re is a coupling introduced to 
control the number of boxes. As shown in [6], the qualitative behavior of the model depends only 
on the sub-exponential factors of the single-box-occupation probability p (q), since a redefinition 

p{q) p' (<?) = e _K ° e^ oq p(q) (3.5) 


just results in 

Z (TV, re) -> Z' (TV, re) = e»° N Z{N,k- re 0 ) . 

For power like sub-exponential weight factors 

p(q) = q~ 0 , ?e N 

it was shown in [7] that the free energy 

F (re) = lim log {Z (TV, re)) 

N^too Tv 


(3.6) 

(3.7) 


(3.8) 


has a singularity at re cr = — log (( (/3)) for [3 £ (1, oo), where ( ( x ) is the Riemann Zeta-function. 
The phase is called fluid for re > re cr and condensed for re < re cr . For /? > 2 the phase transition 
is 1 st order whereas for f3 £ Fzj transition is n th order. The order parameter for the 

transition is the first derivative of the free energy (3.8) with respect to re, which yields the average 
number of boxes divided by the number of balls 


d F(k) (M) 

—-- = lim ——, 

dn n^-oo TV 


(3.9) 


which vanishes in the condensed phase and equals 1 in the fluid phase. 


The relation to the 4-dimensional EDT model is normally established by identifying the triangles 
(or nodes) in the triangulation with boxes and the number of balls in a box with the number of 4- 
simplices that share the corresponding triangle (node). In this way, the coupling K’> of the EDT 
model nicely takes over the role of the re in the balls in boxes model, the average Regge curvature 
becomes in the thermodynamic limit the analogue of the order parameter (3.9), and the /? in (3.7) is 
related [7] to the /3 used in EDT models with a modified measure term [2], Alternatively, one can 
identify the bubbles or baby-universes with the boxes and the number of necks of each bubble with 
the number of balls in the corresponding box [ 8 ], This yields an effective theory for EDT in the 
form of a branched polymer model, in which the bubbles are the vertices and the necks correspond 
to links between the nodes (as in the figures 2 , 8 , but ignoring the different sizes of the nodes). 
While the latter yields just an effective theory, the problem with the former correspondence is, that 
due to geometric constraints, the interplay between the number of 4-simplices per triangle (or per 
node) and the number of triangles (nodes) itself is much more involved than the interplay between 
the balls and boxes in the balls in boxes model. We would therefore like to propose a different 
correspondence in which the numbers of "balls" and "boxes" are less constrained. 

To this end, let us focus on the largest bubble of a triangulation which we will from now on also 
call base-manifold : this largest bubble is made up of elementary building-blocks consisting of 4- 
simplices and minimal necks 9 . Now consider these elementary building-blocks of the base-manifold 

^Instead of thinking of a minimal neck as a kind of worm-hole to a baby-universe, rather think of the baby-universe as the 
blown-up interior of a space-time region that has the boundary of a 4-simplex, i.e. a minimal neck. 
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as boxes and the number of 4-simplices contained in them as balls. An elementary building-block 
that is an ordinary 4-simplex corresponds to a box containing a single ball, whereas an elementary 
building-block consisting of a minimal neck corresponds to a box that contains as many balls as 
there are 4-simplices in the baby-universe branch behind that neck. The 4 dimensional EDT model 
therefore corresponds to a balls in boxes model with N 4 balls, where each box contains at least one 
ball and the number M of boxes can vary from 6 (minimal size for a (combinatorial) simplicial 4- 
sphere) to N 4 (no necks in the triangulation). The canonical EDT partition function could therefore 
be interpreted as the k = 0 case of the more general partition function 


Z (k 2 ,N 4 ,k) 


Ni 


Ni 


Y eKM Z(k 2 ,N 4 ,M) = Y e KM Y CV e 


«2tV2(T) 


M =6 


M =6 


TgT(M,N 4 ) 


(3.10) 


where T (M , N4) is the set of triangulations that possess a largest bubble of size M and consist in 
total of N 4 4-simplices. 




Figure 19: Average (normalized) number (M) /N 4 of elementary space-time building blocks in the largest 
bubble (left) and the corresponding susceptibility (right) as a function of k 2 for systems of total size N 4 « 32k 
(dark blue), N 4 « 48k (dark red) and N 4 « 64k (dark green). 


In terms of (3.10) the average size of the largest bubble and the corresponding susceptibility 
shown in Fig. 19 could for example be expressed as 


<M) (k 2 ,JV 4 ) _ 1 dlnZ (k 2 , N 4 , k) 




N. 


dn 


K =0 


({M — (M)) 2 ^ (k 2 , -^ 4 ) _ 1 d 2 \nZ{K 2 ,N 4 ,K) 


N 4 


N 4 


8k 2 


Kj—0 


(3.11) 

(3.12) 


A histogram for the M /A 4 -distribution at the pseudo-critical point is shown in Fig. 20 and should 
be compared with Fig. 2 of Ref. [6], As can be seen. Fig. 20 looks much more like Fig. 2 of Ref. [6] 
than the N2/N4 -distribution shown in Fig. 7, which would be the corresponding quantity according 
to the old identification: triangles —► boxes, 4-simplices -4 balls. In particular, M/N4 is a nice order 
parameter as it tends to zero in the elongated phase, while N 2 /N,\ remains finite. 
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Figure 20: Histogram for the fraction of 4-simplices contained in the largest bubble of a system of size IV 4 ~ 
32k at K 2 = 1.2577 (left) and = 1.2587 (right). This quantity serves as an approximation for M/N 4 , 
the number of elementary building-blocks of the largest bubble divided by the total number of 4-simplices, for 
which we would have to take into account also the number of necks of the largest bubble, which would lead to 
about 12 % larger values. 


In what follows, we will show that the Z ( k 2 , 7V 4 . M) appearing in (3.10) can be written as 

Ni + l-M / M \ 

Z(K2,N4,M) = Zq(k 2 ,M) e n p(K2,n k ) I ^N 4 ,n 1 -\ -h n M ) (3.13) 

\k =1 / 

where the first factor Zq (re 2 , M) corresponds to the average number of ways a base-manifold con¬ 
sisting of M elementary building blocks (i.e. minimal necks or ordinary 4-simplices) is realized at 
k- 2 , and the second factor, the sum, is the corresponding probability for N 4 4-simplices to fit into the 
M elementary building blocks, with p (k 2 , n) being the probability for a single elementary building 
block, to have volume n. 

To write Zq (k 2 , M) and p (k 2 . n) more explicitly, we need the micro-canonical partition function 
Z\ (N 2: N 4 ) that counts the number of possible triangulations with N 2 triangles, N 4 4-simplices 
and which have a boundary of the form of a minimal neck. For JV 4 > 5 each such triangulation can 
be obtained by removing a 4-simplex from a corresponding triangulation without boundary, that has 
the same number N 2 of triangles but consists of (IV 4 + 1) 4-simplices. Thus Z\ (N 2 , N 4 ) can be 
expressed in terms of the ordinary micro-canonical partition function Z (N 2 . N 4 ) as 10 : 


Zi(N 2 ,N 4 ) 


1/5! , 1V 4 = 1, N 2 = 10 

(N 4 + 1 ) Z (N 2 , N 4 + 1) ,N 4 > 5 


(3.14) 


where 1/5! is the symmetry factor of a 4-simplex and (N 4 + 1) is the number of possibilities to 
remove one 4-simplex from a triangulation of size (JV 4 + 1). 

The corresponding canonical partition function is then: 


Z\ (ft 2 , JV 4 ) = V E e ^(N 2 (T\aT) + ^ 2 N 2 (0T)) 

‘ ^ Gt 

TSTi(AT 4 ) 


J2 Zi(N 2 ,N4)e K ^ N2 ~ 5 \ (3.15) 

n 2 


where T) (N 4 ) is the set of triangulations with a minimal boundary that consist of N 4 4-simplices. 
We can now express Zq (k 2 . M) in terms of (3.15) and (1.6) by noting that the number of ways in 
which M elementary building blocks can be glued together to form a base-manifold, is the same 
as the number of ways to form a triangulation, consisting of M 4-simplices, that does not have any 

10 We neglect complications due to changing symmetry factors CV■ occurring when removing a 4-simplex from triangula¬ 
tions, as the number of symmetric configurations contributing to Z (IV 2 , IV 4 ) is hopefully negligible for large AT 4 . For small 
N 4 , it might be necessary to take the effect of a changing CY into account. 
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neck. We can therefore write: 


Z 0 (« 2 , M) = Z (k 2 , M) — 5! ^ Zi (/c 2 , M — m) Z\ (k 2 , tu) 

m 

= Y (^(A r 2 ,M)-5!^^Z 1 (iV 2 -n2,M-TO) Z 4 (n 2 ,m) ) e K2 JVa , (3.16) 

N 2 \ m n 2 ) 

where the first term within the brackets on the second line corresponds to the number of triangula¬ 
tions consisting of M 4-simplices and N 2 triangles, while the second term, which is a sum over all 
possibilities to form a triangulation of size M by glueing two triangulations, each with a minimal 
boundary, along their boundaries (in [3], this second term was used to measure the entropy exponent 
by baby-universe counting), subtracts the subset of these triangulations that in addition possess at 
least one minimal neck 11 , such that the whole bracket yields the number of triangulations with M 
4-simplices, N 2 triangles and no necks. 

The probability distribution p (re 2 , n) required for the second factor in (3.13) is given by p (re 2 , 7V 4 ) oc 
5!Z lnt (K 2 ,7V 4 ), where 

Z[ nt (n 2 ,N 4 )= Y 7T- e K2Ar2(TV9T) = Y Z i ( N 2 , N 4 )e K2(N2 - 10) (3.17) 

TeTi(Ati) L ' T N 2 

is the canonical partition function for the interior of triangulations at k 2 that consist of iV 4 4- 
simplices and possess the boundary of a 4-simplex. 

The reason for subtracting the whole boundary from the action in (3.17) is, that these terms are 
already taken into account in Z 0 (k 2 , M) and we want to avoid over-counting 12 . 

After having written the EDT partition function in the generalized form (3.10), i.e. in terms of a 
base manifold and its elementary volume elements, which can be excited to form "baby universes", 
some comments are in order: 

1. The terminology "base-manifold" and "elementary building blocks" already suggests that we 
would like to look at the triangulations, observed in EDT simulations, in a slightly non¬ 
standard way. The main reasons for such a re-interpretation are the following: 

• it seems that the base-manifold, with all elementary volume elements in the "ground 
state" (such that they are just ordinary 4-simplices), can be mapped on a corresponding 
Lorentzian or causal triangulation, 

• although the boundaries of the elementary volume elements are always minimal, their 
volume can now change in a discrete manner. This makes EDT to fit a little better into 
the quantum gravity picture provided by spin-foam models. 

2. The altered physical interpretation suggests, that the thermodynamic limit should be taken by 
sending M, the number of elementary building blocks of the base manifold, to infinity instead 
of (just) iV 4 13 . 

3. According to [16], the phase transition is associated with a change of sign in the effective 
curvature. In the crumpled phase, the base-manifold (or "mother universe") has negative cur¬ 
vature: there are two singular vertices, which could be seen as the centres of two hyperbolic 
4-balls (each of them formed by many 4-simplices that are all incident to the same central 

11 Again, corrections due to changing symmetry factors (when gluing two triangulations along their minimal boundary) 
might be necessary in (3.16). 

12 Including the boundary terms as in (3.15) into the action for the elementary building-blocks does not work as now at 
least three boundaries (not just two) meet at each boundary-triangle. The boundary action of an elementary building-block Ti 
is therefore given by S ( dTi ) = k ,2 52 a edT n (A) ’ w ^ ere A runs through all triangles in dTi and n (A) is the number 
of boundaries which contain A. The p(k 2 , N 4 ) would then depend (through the n (A)) on the connectivity of the base 
manifold, which is highly undesirable. 

13 As each elementary building block of the base-manifold contains at least one 4-simplex, M —> 00 implies A 4 —> 00 , 
but the converse is not true. 
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vertex) that are glued along their boundary to form a topological 4-sphere. Without a term 
in the action that prevents the base-manifold from shrinking, it seems to be favourable for 
a triangulation to collapse into a baby-universe tree as soon as the singular vertices disap¬ 
pear. Running simulations at quasi fixed M instead of quasi fixed N 4 (but with re,| sufficiently 
large), or at a non-zero value of the new coupling re, would prevent the triangulation from such 
a collapse. The base-manifold should then survive the disappearance of the singular vertices 
and develop a positive effective curvature itself (instead of generating the positive effective 
curvature by producing many small bubbles), which would give rise to a new phase and a new 
phase transition that could be of higher than 1 st order. 

4. For finite systems, the role played by the new coupling re in (3.10) is related to that of the 
anisotropy factor in CDT as re affects the ratio of the average diameter 14 (~ average time 
needed to pass through) and volume of the elementary building-blocks of the base manifold. 

In a follow-up paper we will try to verify the above assumptions and study the properties of 
(3.10) in more detail. An interesting question is of course whether for some values of re, (3.10) 
yields a 2 nd or higher order transition in re 2 (or the fixed p = N 4 /M, or fixed M version of (3.10)) 
and if, when integrating out the volume fluctuations of the elementary building blocks of the base 
manifold, re and k,\ can be combined to yield a kind of effective cosmological constant, such that one 
recovers the form of the original Euclidean Einstein-Regge action. Alternatively one could interpret 
the additional weight in (3.10) as a measure term. 


4 Conclusion 

Our study confirms the qualitative findings of [5,9]: for K 2 ~ re J 2 " r (IV 4 ) we find for N 4 > 32k a 
clear double peak structure in the IV 2 distribution, which becomes more pronounced with increasing 
system size (and there is no sign that the two peaks will eventually merge again in the thermody¬ 
namic limit). This is characteristic of a weak 1 st order transition. A finite size scaling analysis of the 
4 th order Binder cumulant of the N 2 distribution confirms this further. 

As the phase transition is 1 st order, finite systems should allow for coexisting phases in a neighbour¬ 
hood of the pseudo-critical point. It turned out to be difficult to give a precise criterion to distin¬ 
guish "locally" between the two phases but a candidate could be that bubbles with a neck density 
Pnecks > Packs’ can t> e considered as corresponding to the elongated phase, where p^ ecks is not yet 
known exactly but seems to be around 1/9. Bubbles with p necks < Pnecks wou ld then correspond 
to the crumpled phase. 

Finally we proposed a new correspondence between the EDT and "balls in boxes" models which 
leads to a generalization of the EDT partition function (with an additional parameter re) and a modi¬ 
fied interpretation of triangulations contributing to the EDT partition sum in terms of a largest bubble 
or "mother universe" and its elementary building blocks, which can undergo volume excitations such 
that their interior could also be interpreted as a baby-universe branch. The additional coupling re en¬ 
riches the phase structure of the model which could now possibly contain a 2 nd order phase transition 
line. 

In the Appendix, we propose and motivate a change in the EDT path-integral measure which intro¬ 
duces tunable parameters ( r n ). For appropriate choices of the r n , the order of the phase transition of 
the ordinary EDT model might also change to second order. 
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14 The average diameter of an elementary building-block (i.e. the average time needed to pass through it) consisting of JV 4 
4-simplices and N 2 triangles (counting also the ones in the boundary) can be determined by measuring the average return 
time of a random walk in systems with (N 4 + 1) 4-simplices and N 2 triangles, according to (3.14). 
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A Appendix: Geometric Probabilities and Path-Integral Mea¬ 
sure 

It has recently been suggested [2] that the 1 st order transition of the EDT model could perhaps be 
changed into a 2 nd order transition by a change of the measure in the partition sum (1.5). However, 
none of these attempts has proved successful. Here we motivate and derive a new proposal for a 
measure that could have the desired properties. 


The measure that is normally used is the trivial one for which we have 


P{T) 


e K 2 AT 2 (T)-K 4 Af4(T) 

C~ T 


e K 2 Af2(T)- K4 Af4(T) 


(A.1) 


in the detailed balance equation (2.2). But we could also introduce a measure z ( T ) on the space of 
possible triangulations T such that we would have 


Z (T) Q K -iN 2 {T)-K i N i (T) 

P{ T ) = ~ - -pi - 


z (T) q k 2 n 2 ( t )- k ± n 4 ( t ) 


(A.2) 


To motivate a particular form of the measure z ( T ), assume we are currently in a triangulation 
T which possesses /„ (T) locations where a n-move could be applied. This means that T has 
Y. fn (T) "neighboring" triangulations. Now think of each location in T where a move can be 

n—0 

applied as something similar to a site in an Ising spin system where a spin-flip can occur. But in 
our case, the "spin-flip" consists of the application of a Pachner-move, e.g. a Pachner //-move that 
flips a piece of triangulation, spanned by (5 — n) 4-simplices into one that is spanned by (n + 1) 4- 
simplices but has the same boundary. As the regions where different moves are possible can overlap, 
the flip of one region will in general destroy some of the other regions where flips were possible and 
instead create new ones. We therefore have a fluctuating number of degrees of freedom (see Fig. 21) 
and the system is much more involved than an Ising system. 




K 2 


K 2 


Figure 21: Average total number of degrees of freedom (i.e. number of possible moves) per volume (left) 
and the corresponding susceptibility (right) as functions of k ,2 for different system sizes N± = 2,..., 64k. 
The red dot on each graph marks the corresponding pseudo-critical point/value (i.e. where k ,2 = K^ cr (IV 4 ), 
corresponding to the peak in the A^-susceptibility). The simulation was performed with the update scheme 
described in Sec. 2.2. 


Nevertheless, one could argue that, as long as the couplings K 2 and K 4 are turned off, all possible 
moves in a triangulation T should be considered as equally likely, just as in the Ising case. More 
generally: one could assign different probabilities r n to different move types n £ {0,..., 4}, as long 
as r n = V 4 - n . For example, one could choose r n proportional to the local volume of a n-simplex 15 
that allows for a n-move (the local volume is the same for n and (4 — n)-simplices which allow for 
a move). 

IS As mentioned earlier: the local volume of a n-simplex is the volume of all points which are closer to this n-simplex than 
to any other n-simplex. 
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Having such fixed probabilities for the different moves, implies that triangulations with different 
numbers of locations where moves could be applied, are not equally likely. The corresponding 
probability weight for a triangulation T can be derived from the balance equation 


= E E 

ra=OT'Gnbr„(T) 


z (T) r n 
T,rmfm{T) 


4 


E E 

n =0 T> G nbr„ (T) 


z (TQ r 4 _ n 

E r m fm(T’y 


(A.3) 


where nbr„ (T) is the set of all triangulations that can be obtained from T by a n-move. The detailed 
balance equation corresponding to (A.3) reads 

z (T) r n _ z (TQ r 4 - n 

E r n f n (T) E r n f n (T') 

n—0 n—0 

which is obviously satisfied if we have 

4 

z (T) cx E r « fn ( T) . 

n=0 

Such measures are particularly simple to implement: by choosing move candidates according to the 
selection probability, 

P { n sl) ( T) = - 4 -^-, (A.6) 

E r m f m (T) 

771 = 0 

the measure term drops out of the detailed balance equation for the (reduced) transition probabilities 

Pn (T): 

p (T) ( T) p n (T) = p (T') (TO p 4 -n (TO , (A.7) 

where p (T) is given by (A.2), such that there is no need to determine the possible moves of the 
candidate configuration. 


(A.4) 


(A.5) 
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